T * H • E 

OHIO 

SME 


UNIVERSITY 


/As ~ / t- - C/Z~ 

/ 3-^d 

Alt 


Material Parameter Measurements at High 

Temperatures 


by 


A. Dominek, A. Park and L. Peters, Jr. 


The Ohio State University 

ElectroScience Laboratory 

Department of Electrical Engineering 
Columbus, Ohio 43212 


(BASA-CB-1827C7) HAT EBIAL E AE ABE1EB N88-2C3S7 

F.EASUBEBEN55 A1 EJGH 1EHPEBA1CBE5 {Ohio 
itate Dniv.) 2A f CSCL 07D 

Onclas 

G3/25 0135352 


Technical Report No. 719300-2 
Grant No. NAG3-784 
March 1988 


National Aeronautics and Space Administration 
Lewis Research Center 
21000 Brookpark Rd. 

Cleveland, OH 44135 




NOTICES 


When Government drawings, specifications, or other data are 
used for any purpose other than in connection with a definitely 
related Government procurement operation, the United States 
Government thereby incurs no responsibility nor any obligation 
whatsoever, and the fact that the Government may have formulated, 
furnished, or in any way supplied the said drawings, specifications, 
or other data, is not to be regarded by implication or otherwise as 
in any manner licensing the holder or any other person or corporation, 
or conveying any rights or permission to manufacture, use, or sell 
any patented invention that may in any way be related thereto. 



Contents 


List of Figures iv 

1 Introduction 1 

2 Alternate Measurement Fixtures 3 

3 Parameter Extraction 7 

4 Canonical Sample Geometries 10 

5 Example Calculations 14 

6 Conclusions 19 


PRECEDING PAGE BLANK NOT FILMED 


iii 

fig _ i #" .wwf nwwwny now* 



List of Figures 

2.1 Illustration of the parallel plate test fixture 4 

2.2 Illustration of the ogival shaped ground plane test fixture. . 6 

4.1 Geometry for a coated cylinder 11 

4.2 Geometry for a coated hemisphere 13 

5.1 Echo widths for circular cylinders. Perfectly conducting - 

solid, homogeneous dielectric - short dashed and coated - 
long dashed. Perfectly conducting and coated cylinder re- 
sults are very similar for these parameters 15 

5.2 Echo areas for hemispheres on an infinite ground plane. Per- 

fectly conducting - solid, homogeneous dielectric - short dashed 
and coated - long dashed 16 

5.3 Sphere diameter sensitivity example to estimate the relative 
permittivity from calculated values based upon a diameter 
of 3” and relative permittivity of 2.6. Assumed diameter: 3” 

- Solid, 2.97” - dot-dashed and 3.03” - short dashed 17 

5.4 Numerically determined relative permittivity from measured 

backscattered data for a 3” diameter sphere 18 


IV 



Chapter 1 
Introduction 


Conventional techniques for the evaluation of the constitutive material pa- 
rameters involves material samples filling an interior region of a coaxial or 
waveguide line. Care must be exercised in acquiring the sample to insure 
a firm fit in the line to obtain accurate reflection and transmission mea- 
surements. The use of such a line fixture to perform these measurements 
becomes undesirable at high temperatures (greater than 800° F) due to 
mechanical and electrical considerations of the line. Alternate approaches 
are suggested here to eliminate the need of closed system measurements to 
acquire material parameters. The upper temperature restriction in these 
approaches would be limited by the fixture durability and the means to gen- 
erate the desired temperatures within these fixtures. The alternate system 
involves the measurement of scattered fields from canonical sample geome- 
tries using parallel plate or ground plane fixtures. The material parameters 
can then be found from measurements and exact numerical expressions for 
the scattered fields from these sample canonical geometries. Similar free 
space measurements have been performed for material coated plates based 
upon a reflection model [2]. 

The swept frequency measurements acquired with either of these fix- 
tures can readily be performed with the HP 851 OB Network Analyzer. The 
software time gating of the HP 8510B would be implemented to eliminate 
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error terms from the edge scattering of the finite fixture size. 
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Chapter 2 

Alternate Measurement 
Fixtures 


There are two possible fixtures that would readily allow the necessary scat- 
tered field measurements from heated sample geometries. They are a paral- 
lel plate and ground plane geometries shown in Figures 2.1 and 2.2, respec- 
tively. The canonical sample geometries used in these fixtures are circular 
cylinders for the parallel plate fixture and hemispheres for the ground plane 
fixture. These geometries are the most natural because their exact scat- 
tered fields are readily computed from eigenfunction solutions for either 
homogeneous material geometries or material coated perfectly conducting 
geometries. 

The parallel plate geometry is ideal for low frequency measurements 
where the wavelength of the incident frequency is greater than twice the 
distance between the parallel sides of the fixture. This restriction is so 
that only one dominant mode exists in the fixture. It is suggested that the 
material sample be coated on a hollow perfectly conducting cylinder to allow 
internal heating by an electric or gas source. This internal heating approach 
would eliminate possible mechanical and electrical difficulties caused by a 
heat source in between the parallel plates. For the heat source between 
the parallel plates, the sample would have to be rotated to insure uniform 
heating. The difficulty arises from the crack that would be formed to allow 
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the rotation. The scattering from the crack would have to be accounted 
for in the analytic model which would be difficult. The rotation of the 
sample would also have to be very uniform along its axis. Another benefit of 
internal sample heating is that no consideration is required for the scattering 
from the heat source. Permitting the heat source between the parallel 
plates would introduce another timing constraint for the separation distance 
between the sample and heat source. The expansion of the sample due to 
the heat is of no great concern as long as it is known for the parameter 
determination. It is necessary however to maintain tight contact between 
the sample and the sides of the parallel plates. 


ACCESS PORT FOR HEAT SOURCE 



TOP VIEW 


Figure 2.1: Illustration of the parallel plate test fixture. 

Many of the 6ame concepts for the parallel plate fixture are applicable 
to the ground plane fixture. Again it would be desirable to heat the sample 
internally to exclude scattering from the heat source and the crack in the 
ground plane which would permit rotation. It would be very desirable to 
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shape the ground plane to have an ogival contour to control the scattered 
field from its edges. The heating hardware would have to be shielded by 
extending a curved metal surface from the ground plane surface. 

It is important that the distance between the test sample and the source 
be large enough to produce a planar illuminating field for the test sample. 
This is for a simplification of the calculation of the exact scattered fields, 
though it is not absolutely necessary. This minimal distance, called the 
range distance R, is given by 


R = 



where £ is the maximum length of the test sample geometry and A is the 
wavelength of the incident field [1]. The maximum phase variation of the 
incident wavefront at this distance is 22.5°. 

The range distance is not the dominating factor for the parallel plate 
fixture if the upper frequency is set at 18 GHz. The distance between the 
parallel plates for only the dominant mode to exist is approximately .325 
inches. The main test bed of the parallel plate fixture should be at least 
1.5 feet square with the sample in the middle. The fixture requires a feed 
region that would consist of a tapered parallel plate geometry expanding 
from a vertex feed to the width of the main fixture. The length of the 
feed should be approximately 1.5 feet long. This fixture would be very 
desirable in terms of providing a uniform, high temperature environment 
to the sample. The major disadvantage for this fixture is that the material 
sample length could be too short for an adequate field interaction to occur. 

The minimal range distance for the ground plane fixture is 5 meters if 
an upper frequency of 18 GHz is desired. A tip to tip length for the ground 
plane should be in the vicinity of 6 feet with a ground plane width of 1 
foot. 
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SIDE VIEW 

SHAPED GROUND PLANE 
FOR MINIMAL EDGE RETURN 
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TO SHIELD 
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TOP VIEW 


Figure 2.2: Illustration of the ogival shaped ground plane test fixture. 
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Chapter 3 

Parameter Extraction 


Extracting the constitutive parameters from measured swept frequency 
data requires two items. One item is calibrated data and the other is an 
exact solution for the measured scattered field. Calibrated data is obtained 
from three measurements, these being test sample data, background data 
(no test sample present) and data from a known reference. This data is 
then used in the following formula to obtain calibrated data 


CSD = 


/ MSP - MBD \ 
VMRD - MBD/ 


CRD 


(3.1) 


where CSD is calibrated sample data, MSD is measured sample data, MBD 
is measured background data, MRD measured reference data and CRD is 
calculated (exact) reference data. The vector subtraction using the mea- 
sured background is not necessary if stray scattered signals are negligible 
or the HP 8510B Network Analyzer time gating window is small enough to 
include only the desired signal. 


Once the measured data has been calibrated, the Newton- Rap hson tech- 
nique can be applied to numerically search for the proper constitutive pa- 
rameters of the material, t and fi. Since these parameters are complex, two 
independent measurements have to be performed. In cases where the pa- 
rameters are not highly sensitive to frequency, two measurements at slightly 
different frequencies will suffice or an average of values around the frequency 
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of interest will also work [2]. Such an approach has to be taken for monos- 
tatic (transmitting and receiving horns are at the same location) scattering 
measurements using the parallel plate fixture. Independent measurements 
at the same frequency for this fixture can be obtained if bistatic measure- 
ments are taken. The ground plane fixture has two extra features available 
for independent measurements which are the polarization of the incident 
field and the elevation angle of the incident field direction from the ground 
plane surface. 

Another potential means to obtain independent scattered field relation- 
ships from one measurement is to isolate the individual scattering mech- 
anisms that occur in the total scattered field. The scattered field for a 
homogeneous dielectric sphere has several wave interactions. The two most 
apparent ones are the reflection from the front surface and the internal re- 
flection from the rear surface. These two reflections are dependent upon 
the constitutive parameters in different functional relationships. Effectively 
two different measurements can be obtained from one measurement if these 
reflections are separated in time enough so they can be isolated from the 
total scattered field by the time gating software in the HP 8510B. 

The Newton-Raphson technique iteratively solves for the unknown pa- 
rameters by forcing two convenient functionals to zero. Defining the func- 
tionals, 1 * 0 , 1 , to be 


F 0 = Mq — Co — 0 


(3.2) 


and 


Fi = Mi - Ci = 0 (3.3) 

where the Afi i2 and Ci >2 are the independent measured and calculated re- 
sponses for the sample geometry. The calculated responses are dependent 


8 



upon the parameters e and /i which are iteratively varied until some mini- 
mum error criterion has been achieved. The unknown parameters are iter- 
ated with the following process from some initial guess 


£i + l — £{ + A 


(3.4) 


and 


Vi+i = Pi + A fii 


(3.5) 


where 


Ae = 


F,2Esi l 

ri 8c 


jp 8F i 

^0 8 e 


(3.6) 


and 


A n ~ 


jp 8Fi jp 8Fn 

*0 - r 1-JT7 


fi/J 


and 


8Fn 

8Fn 

8c 

8n 

8Fi 

8 Ft 

8c 

8n 


(3.7) 


(3.8) 


9 


Chapter 4 

Canonical Sample Geometries 


The exact scattered fields for the canonical sample geometries are readily 
calculated by eigenfunction expressions. The geometry for a coated per- 
fectly conducting circular cylinder is shown in Figure 4.1. The scattered 
far field for plane wave illumination is given by [3] 

/ O p ~jkp oo 

; e ’' 5 7^£™ ( 2-* ) (4 ' 1) 

where a n are ratios of cylindrical Bessel functions, J n , N n and and 
given by 

Jn(M) ~ Z n J n (kb) 

en H i 2) \kh) - Z n Hn ] (kb) 

and 


r . jjMmw - (M) ,,,, 

" V ’ Mk,a)N4k 1 i) - N„(k ia )J„(k,k) { l 

where rj r = is the relative wave impedance of the material and fcj = 
y/ € T fJL r k and k is the free space wave number. The parameter e n is equal to 
unity for n = 0 and 2 for n ^ 0. 

These above expressions reduce to simpler forms for special cases. The 
quantity Y n is equal to 
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-rW) 

Vr Jnikxb) 

for a homogeneous cylinder. Additionally, for the perfectly conducting 
cylinder, Y n — oo. 



Figure 4.1: Geometry for a coated cylinder. 

The geometry of a coated perfectly conducting hemisphere over an infi- 
nite ground plane is shown in Figure 4.2. The scattered far field for plane 
wave illumination can be calculated by using image theory which results in 
the sum of the backscattered and bistatic scattered fields for a sphere in 
free space. The scattered field for a hemisphere on a ground plane can be 
written in the following form for the two principal polarizations 


e, = = o) + = 2#„) 


(4.4) 


and 


E+ = E ,phere (6 = 0) - E‘/ here (8 = 20 o ) (4.5) 

where E‘ e p £ ere is the scattered field for a sphere in free space at a bistatic 
angle of 9 a . 
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The scattered field expressions for a sphere in free space are given by 


[ 3 ] 


E; phere (8) 


171 oo 

-ir e ' ikr E 

kr n^l 


b n sin 0P „( cos 6) 


P^cosfl) 
n sin 0 


( 4 . 6 ) 


and 


E?"V) = 


E 



Pn( COS 0 ) 

sin 9 


c n sin flP^cos 6) 


( 4 . 7 ) 


where P^cosfl) and P 1/ (cos5) are Legendre polynomials and 



B n {A n F\ - t h X n Fl) 

A n {B n FX-VrX n FX) 


( 4 . 8 ) 


and 


Cfl 


. (X„F„ 3 -,,X„F*) 
\X.F* - n,B„FZY 


( 4 . 9 ) 


where r] r = is the relative wave impedance of the material. The 
symbols 6° , c° , A n , B n ,X n and F^’ 2,3 ’ 4 are ratios of spherical Bessel functions 
j n and h n given by 


2n + 1 [kbj n (kb)}' 
n " n ( n + 1) ^kbhn\kb)^'' 

(4.10) 

0 2n + 1 jnW 

n(n + l) *!?>(*»)’ 

(4.11) 

, |Wj„(B)f 

kbj n (kb) 

(4.12) 

[»*!?)(»)]' 
n kbhL 2) (kb ) 

(4.13) 
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Min(fcl&)r 

kibj n {k x b) 


[hajniha)]' 

jn(kib) |fc 1 a/4 l 2 ^(fcia)j 

[kiajn(kia)]' [ k i b ^ 2 *(&!&)] 

[kibj n (kib)] 

*,(*,.) \k,ihP(k,i)\' 

[Mjn(M)F h£\k t a) 


and 


jn(fcia) hWfab) 
j n (kib) h^i^a) 


where ki = y/e r (i r k and k is the free space wave number. 

These above expressions also reduce to simpler forms for special 
The F*' 2,3,4 quantities are equal to unity for a homogeneous sphere, 
tionally, for the perfectly conducting sphere, b n = 6° and c n = c°. 



Of ) 



Figure 4.2: Geometry for a coated hemisphere. 


(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 
cases. 

Addi- 
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Chapter 5 

Example Calculations 


The examples that follow have been chosen for simplicity and ease of com- 
putation to illustrate the concepts involved. The thin material coatings are 
lossless with fi r = 1.0. Obviously the introduction of thicker, lossy coatings 
with fi r ^ 1.0 will increase the influence of the material on the scattered 
fields from the conducting surfaces. These solutions should be carried to 
a conclusion in further studies in a manner similar to the work of Yu and 
Peters [4]. 


The backscattered fields of the canonical geometries are calculated for 
comparison purposes between the perfectly conducting case, the homoge- 
neous dielectric case and the dielectric coated-perfectly conducting case. 
The calculated fields are expressed as echo widths for the circular cylinder 
geometries, given by 


<tw — lim 

p—* OO 


E' 2 

OF 


(5.1) 


and echo areas for the hemispherical geometries, given by 

, r E ‘ '* 

<7 a = 4n lim — 

r-oo E l 


(5.2) 


where E *•*’ are the scattered and incident electric fields. Figure 5.1 shows 
the echo widths for a 2” diameter perfectly conducting circular cylinder, 2” 
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diameter dielectric circular cylinder with e T = 2.6 and a coated 2” diameter 
perfectly conducting circular cylinder with a coating thickness of .04” and 
e r = 2.6 calculated between 2 and 18 GHz. In this case, the coating for its 
electrical thickness and constitutive parameters has negligible influence on 
the scattered fields in this frequency range. Its scattered fields for the coated 
case are essentially the same as the perfectly conducting case. Figure 5.2 
shows the echo widths for a 3” diameter perfectly conducting hemisphere, 
3” diameter dielectric hemisphere with e r = 2.6 and a coated 3” diameter 
perfectly conducting hemisphere with a coating thickness of .04” and e r = 
2.6 calculated between 2 and 18 GHz. In all calculations, fi r = 1. 



FREQUENCY IN GHZ 


Figure 5.1: Echo widths for circular cylinders. Perfectly conducting - solid, 
homogeneous dielectric - short dashed and coated - long dashed. Perfectly 
conducting and coated cylinder results are very similar for these parameters. 


These figures readily illustrate the minor influence the material coating 
has on the scattered fields when the material parameters are not signifi- 
cantly different from free space. The influence of the coating becomes more 
pronounced with thicker electrical coatings. The homogeneous dielectric 
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Figure 5.2: Echo areas for hemispheres on an infinite ground plane. Per- 
fectly conducting - solid, homogeneous dielectric - short dashed and coated 
- long dashed. 

cases do not have this limitation due to the absence of the inner conduct- 
ing core. 

A sensitivity study was performed on a numerical technique to extract 
the relative permittivity, e r from calculated backscattered fields from a 3” 
diameter homogeneous dielectric sphere with e r = 2.6. A functional similar 
to that described earlier was used except calculated data was used instead 
of measured data. Figure 5.3 shows the calculated relative permittivity 
assuming the sphere diameter to be 2.97”, 3.0” and 3.03”. A 1% change 
in the assumed sphere diameter created approximately a 2% error in the 
calculated relative permittivity. 

The final example is the determination of the relative permittivity of 
a relatively lossless 3” dielectric sphere from backscattered measurements. 
Figure 5.4 shows the results between 2 and 10 GHz. A value of e r = 2.6 
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Figure 5.3: Sphere diameter sensitivity example to estimate the relative 
permittivity from calculated values based upon a diameter of 3” and relative 
permittivity of 2.6. Assumed diameter: 3” - Solid, 2.97” - dot-dashed and 
3.03” - short dashed. 

was obtained which is reasonable for the material. 
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Figure 5.4: Numerically determined relative permittivity from measured 
backscattered data for a 3” diameter sphere. 
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Chapter 6 
Conclusions 


The constitutive parameters can be obtained in principle from scattered 
field measurements of a material body. The ground plane fixture required 
for such a technique does not involve the temperature related difficulties 
which the traditional technique for material parameter measurements has. 
Of the two scatterers examined, the hemispherical geometry is more de- 
sirable in terms of fixture design and electrical performance. The most 
reliable parameter determination from measured scattered fields of a mate- 
rial coated, perfectly conducting scatterer is obtained with electrically thick 
coatings. 
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